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On the dust abundance gradients in late-type galaxies: 

I. Effects of destruction and growth of dust in the interstellar 

medium 
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ABSTRACT 

We present basic theoretical constraints on the effects of destruction by supernovae (SNe) and 
growth of dust grains in the interstellar medium (ISM) on the radial distribution of dust in 
late-type galaxies. The radial gradient of the dust-to-metals ratio is shown to be essentially 
flat (zero) if interstellar dust is not destroyed by SN shock waves and all dust is produced 
in stars. If there is net dust destruction by SN shock waves, the dust-to-metals gradient is 
flatter than or equal to the metallicity gradient (assuming the gradients have the same sign). 
Similarly, if there is net dust growth in the ISM, then the dust-to-metals gradient is steeper than 
or equal to the metallicity gradient. The latter result implies that if dust gradients are steeper 
than metallicity gradients, i.e., the dust-to-metals gradients are not flat, then it is unlikely dust 
destruction by SN shock waves is an efficient process, while dust growth must be a significant 
mechanism for dust production. Moreover, we conclude that dust-to-metals gradients can 
be used as a diagnostic for interstellar dust growth in galaxy discs, where a negative slope 
indicates dust growth. 

Key words: Galaxies: evolution, ISM; ISM: clouds, dust, extinction, evolution, supernova 
remnants; 



1 INTRODUCTION 

The lifetime of dust grains in the interstellar medium (ISM) is a 
critical parameter for the evolution of the dust component in a 
galaxy. Shock-waves originating from supernovae (SNe) arguably 
contain enough energy to destroy (or at least shatter) dust grains as 
these waves propagate through the ISM. The time scale for such 
dust destruction depends on several physical conditions, where 
the supernova rate (SNR) and efficie ncy of dust d estruction in a 
SN-shock are the most important jMcKeej fl989l : iDrainel 1990). 
Shock destruction of d ust grains has been considered quite e f ficient 



late-type galaxies seems to work nicely with out significant net de- 
struction of dust dlnouell2003l ; lHirashita|[l999h . In fact, no or little 
net dust destruction makes it easier to explain the dust-to-gas ra- 
tios, since stellar dust production is likely not sufficient for neither 
the Milky Way, nor any of the other late-type loca l group galaxies 
jHirashitalll999l : IZhukovska. Gail & Trielofl1l2008h . 

Observation al constraints imply du st production in SNe is 
rather inefficient dKotak etal.ll200al2009h . which suggest the high 
dust masses detected in some, relatively old, SN remnants (see. 



in many studies, e.g., Ijones. Tielens & Hollenbachl Jl996h : Ijonesl 120091 : iGall. Hiorth & Andersen ] 1201 lb could be the result of sub- 



Morgan & Edmunds! 120031; iMorgan et al J 120031 : iDunne et al.l 



d2004h :ls erra Dfaz-Cano & Jonesl ( 12008b . suggest a grain lifetime 
of a few times 10 s yr for many d ust species, but note that a recent 
re-evaluation of dust lifetimes by I Jones & Nuthl J201 1) showed that 
the lifetimes of silicate grains may be comparable to the injection 
time scale of such grains. 

While this high dust-destruction efficiency seems consistent 
with the Milky Way (solar neighbourhood), it has been shown by 
several authors that very efficien t dust destruction is unlikely in 
high-z objects toweketal1l2007l ; IGall. Andersen & Hiorth1l201ll : 
iMattssonlHoITI) . It may of course be that dust destruction by SNe 
is less efficient in high-z galaxies, but also modelling of nearby 
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sequent dust growth, and/or heating of pre-existing dust, rather 
than dust production in the actual SN. However, since SN rem- 
nants with large dust masses typically contain vast amounts of cold 
dust, some degree of growth appear to be necessary even if there 
is a component of heated swept-up dust. This picture is consistent 
with theoretical results which suggest 90% of the dust produced in 
SNe is destroyed by the rev erse shock before it reaches the ISM 
dBianchi & Schne ider 2007). Hence AGB stars could be the most 
important source of stellar dust as a significant fraction of the met- 
als expelled by these sta r s is expected to be in the form of dust 
jEdmunds & Ealeslll998l : iFerrarotti & Gaill I2006T) . which is sup- 
ported by observational detec tions of dust (see, e.g., the recent re- 
sults by Ramsted t et al.l201lb . It should be noted, however, that this 
picture may need to be revised due to the recent discovery of a large 
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amou nt of cold dust associated with SN 1987A jMatsuura et al.l 

l2Qiih . 

In mode ls of dust evolution for the solar neigh bourhood by 

iDweld 1 19981) and IZhukovska. Gail & Trielofll fc008l) . the limited 
dust production in stars and possible dust destruction by SN-shock 
waves are more than well compensated by an efficient dust growth 
in the ISM, which is supported by observations indicating the ex- 
istence of large, micromet er-sized dust grains in dense molecular 
clouds jPagani et alj|201ol) . T here is further evidence from dust- 
evolution modelling (see, e.g. lMattssonll20lH ; IPipino et alj|201 ll ; 
IValiante et al.l 1201 1[) along with some observational constraints 
(see, e.g. Michalowski et alj|2010h which suggests a need for sig- 
nificant dust growth in the early Universe as well. Dust growth ap- 
pears to domi nate over dust destruction also in the local, present- 
day Universe jHirashitjl 19991: Ilnouell2003l ; Iffirashita & Kucll201ll ; 
lAsano et al.l2012h . It is difficult to separate one scenario where dust 
growth is totally dominating over dust destruction from another 
where there is very little dust destruction and less dust growth. But 
whether there is net growth/destruction it should have observable 
consequences, however. 

We propose here a diagnostic tool for determining whether 
there is net dust growth or net dust destruction in the ISM of a late- 
type galaxy for which dust-to-gas as well as metallicity gradients 
can be derived. As we will show in this paper, the change (gradient) 
of the dust-to-metals ratio along a galactic disc is closely connected 
to g rowth and destruction of dust in the ISM. In an associated pa- 
per dMattsson & Andersen! l201li hereafter cited as Paper II), we 
investigate the implications of observed dust-to-m etals profiles in 
a sele ction of late-type galaxies from the SINGS dKennicutt et al.l 
120031) sample. 



2 BASIC EQUATIONS 

In order to obtain analytical solutions and be able to manipulate 
the basic equations of the dust-enrichment problem in such a way 
that we can derive some basic constraints, we use the instantaneous 
recycling approximation (IRA, which essentially means all stars 
are assumed to have negligible lifetimes with respect to the over - 
all time scale for the build-up of metals and dust, see IPagelll 19971) 
throughout this paper. No delayed element production due to stellar 
lifetimes is considered. 

For convenience, we also define the dust destruction rate rela- 
tive to the growth rate of the stellar component S s as 



D(r,t) = l,i SM {r,t)\—j 



(1) 



where %sm is the dust destruction rate due to SNe and the vari- 
ables r, t are the galactocentric distance and time/age respectively 
(a notation that we will assume is understood in the following). 
Similarly, we also define 



(2) 



where £ gr is the rate of grain growth (in mass units) in the ISM. 

Assuming a 'closed box', where dust destruction in the ISM 
is from SN-shocks, the equations for the metallicity Z and the dust- 
to-gas ratio Z d becomes 

dZ dZ s cffi g 

-yz-r, (3) 



dZi dL. cTL. 
, — - = y d — ^ + Z A (r, t) \G(r, t) - D(r, t)] 
at at 



dt ' 



(4) 



where E g is the gas density, X d is the dust density, and the yield y, 
is defined as 



i = — I p,(m)m(f>(m)dm, 



(5) 



for both metals (i = Z) and stellar dust (i = d). In equation l[5} 
above, /?, is the fraction of the initial mass m of a star ejected in 
the form of newly produced metals or dust, a is the stellar lock-up 
fraction (i.e., the fraction of the baryon mass being locked up in 
long lived stars) and <p(m) is the mass-normalised IMF, with ni\ , 
m up being the lower and upper mass cuts, respectively. Combining 
equations {3]l and I0, we have 



dZi = y A +Z A [G{r,t)-D(r,t)] 
dZ y z 



(6) 



which thus have no explicit dependence on the gas mass density E e 
or the stellar mass density E s , although G, D and Z d of course may 
have implicit dependencies on the amount of gas and stars being 
present in a certain galactic environment. 



3 CONSTRAINTS ON DUST-TO-METALS GRADIENTS 

We will now prove some basic properties of dust-to-metals (£) 
gradients relative to the metallicity (Z) gradient. For 'logarithmic' 
dust-to-metals and metallicity gradients we use the following nota- 
tions, 



31nZ d 
= dr = 

dWZJZ) 
dr 



Z d dr ' 

31nZ d 
~ dr 



<91nZ _ 1 dZ 
dr ~ Z~dr' 

1 dZ A 1 dZ 

dr Z A dr Z dr 



dlnZ 



(7) 
(8) 



which are used since they both have the same unit ([length] -1 ). The 
two gradients A f and A z can be regarded as coupled through a func- 
tion / which may be seen as a function of a number of physical 
parameters, but in general we may say it is a function of time t and 
radial position (galactocentric distance) r along the galaxy. Hence, 
we consider a relation of the form A f (r, t) = f(r, t) A z (r, t). In the 
following we will implicitly assume all quantities except y d , yz are 
functions of r and /. We will also refer to the case of a zero deriva- 
tive with respect to r as a 'flat' gradient, which of course could be 
seen as the case of no gradient. However, we prefer to describe the 
gradients as being either positive, flat or negative, where 'negative' 
refers to a gradient (derivative) which decreases with galactocentric 
radius and vice versa for 'positve' gradients. Below we also use the 
sign function sgn(x) = x/\x\ to denote the sign of A f and A z . 

THEOREM. For a closed-box model, without any pre-enrichment, 
and where the IRA and constant yields y z , y d have been adopted, 
the following always hold: 

(i) A flat (no slope) dust-to-metals gradient can only be obtained 
if there is neither net growth, nor any net destruction of dust in the 
ISM (G = D) or if the metallicity gradient is flat. 

(ii) If the dust-to-metals and metallicity gradients have the same 
sign, there has to be net growth (G > D) of dust in the ISM. 

(iii) If the dust-to-metals and metallicity gradients have opposite 
signs, there has to be net destruction (G < D) of dust in the ISM. 



dt 



■-yz- 



dt 



dt 
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Proof. From the basic equations of dust evolution (see section [2} 
one finds 



dZ y z 



(G- D) 1 

yz 



By use of the chain rulfl we get 



d£ 
dr 



yz Z \ yz 



dZ 
~dr~' 



(9) 



(10) 



which in terms of Kr and A z , can be written as 



Vd 1 Z 

— - + — (G-D)- 1 

yz i yz 



The function / (see definition above) is then generally expressed 

(12) 



ya 1 z 

f=—- + —{G-D)-l. 

yz $ yz 



(i) If A z = 0, then obviously A f = as a consequence of Equa- 
tion 1 11 It . In case there is neither net growth, nor any net destruction 
of dust in the ISM (G = £>), we have 



yz <T 



(13) 



Equation ® gives 
dZ A _ y^ 
dZ y z ' 

and again by the chain rule, 

/ dZ A \ _ y d dZ 
I dr j G=D y z dr ' 

Integrating equation d 1 5b . together with the natural initial condi- 
tions Z(r, 0) = Zd(r, 0) = (no pre-enrichment), one obtains 

£ = yjyz, or 
yd 1 



(14) 



(15) 



yz £ 



= l. 



(16) 



Hence, according to Equation d 1 3b . we must have Af = 0, since 
/ = 0, which proves part (i). 

(ii) First, we note that if sgn(A <r ) = sgn(A z ), then / > 0. In case 
G > D, Equation ® gives 



(dZi) 
\ dr 



dZ A 
dr 



(17) 



/ G>D V " IG=D 

Then, by Equation ([8]l and the fact that / = if G = D, we have 
Af,G>£> > A f .o =0 = 0, which implies / > 0. In case G < D, Equa- 
tion ® gives 



(dZi) 
{ dr ' 



< 



dZ A 
dr 



(18) 



lG<D \ IG=D 

Again, using Equation l[8]l and / = if G = D, we have 
Af,o<D < A f _G =0 = 0, which implies / < 0. Hence, 

' If £ is a function of r and t with continuous first partial derivatives, and if 
r and t can be regarded as diff erentiable functions of Z, then 

d( _ d( dr d( dt 

dZ~ drdZ + ~didZ' 

At a specific time t = to we can thus write 

d{ _ d( dr 
dZ dr dZ 



sgn(A f ) = sgn(A z ) is possible if (and only if) G > D, which 
proves part (ii). 

(iii) In this case, if sgn(A^) t sgn(A z ), then / < 0. If G > D, 
Equation ((6) gives 



dZi 
dr 



dZ A 
dr 



(19) 



Analogous to case (ii) we have A ( G>D > A f G=D = 0, which implies 
/ > 0. In case G < D, Equation ® gives 



'II- |^| < | ^ 



dZ^ 
dr 



(20) 



G<D \ IG=D 

Thus, we have A f G<D < A f G=D = 0, which implies / < and 
therfore sgn(A f ) t sgn(A z ) is possible if (and only if) G < D, 
which proves part (iii). 



4 SIMPLE MODELS OF DUST GROWTH AND DUST 
DESTRUCTION 

4.1 Dust growth in the ISM 

The most likely dominant type of 'secondary' dust production is 
that by accretion of atoms (or small molecules) onto pre-existing 
interstellar dust grains. Dust grains can in principle also grow by 
coagulation, but this process will not affect the total dust mass very 
much since it is mostly smaller dust grains being joined together 
into larger grains. Hence, we will here only discuss dust growth by 
accretion. 

We define the rate per unit volume at which the number of 
atoms 7V A in dust grains grows by accr etion of met als onto these 
dust grains in a similar way as (see, e.g. |Pweklll998h 

dN A 



dt 



■ / s na « z « gl (v g ), 



(21) 



where n z and n gr are the total atomic metals and dust-grain number 
densities in the ISM, respectively, a is the typical grain radius and / s 
is the sticking coefficient (i.e., the probability that an atom will stick 
to the grain). (v g ) is the mean thermal speed of the gas particles 
(including metals), which is defined as 



<v 



>>-r 



vf(v)dv ■■ 



MT 



(22) 



where /(v) is the Maxwell distribution, k is the Boltzmann constant, 
T is the kinetic temperature of the gas and m A is the atomic weight 
of the gas particles. In terms of surface densities in the molecular 
gas clouds where the dust may grow, we can write 



d^i 
dt 



(m„)d c 



(23) 



where E z is the surface density of free (atomic) metals, (m gr ) is 
the mean mass of the dust grains in the ISM and d c is the size of 
the molecular cloud in which the dust is growing. The timescale of 
grain growth can then be expressed as 



where 



(•-I)" 



(mgr) d c 



(Wgr) dc 

fs na 2 I.z(v„) f s 7M 2 ZE mol (v g > ' 



(24) 



(25) 
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in which S mo] is the surface density of molecular gas, and Z the 
metallicity. 

For simplicity we will assume Y. mol « E H , , since most of the 
gas in the molecular gas clouds is in the form of molecular hydro- 
gen. We also assume Eh 2 traces the star-formation rate, i.e., 



1 dZ s 

■ ^a 2 = — -T-, 

a at 



(26) 



as ind ic ated by several o bserva ti onal studies ( e.g., Rownd & Young 



19991 : IWong & BlitJ |200A H 
2008; Bigiel et alj |201l[ 



gieletal.1 120081: iLerov et al 



Feldmann, Gnedi n & Kravtsov 



20111 : ISchruba et al.1 120111) 
supported by theory and 



Such a 
recent 



relation 
numerical 



ments 



(see, 



is also 
experi- 



Krumholz, McKee & Tumlinsor] [2009; 



iKrumholz. Lerov & McKeej " 201 1 ). Moreover, the 



mal speed (v g ) is roughly constant in the 



mean 
considered 



ther- 
ISM 



environment and the typical grain radius does not vary much. 
Hence, the timescale To is essentially just a simple function of the 
metallicity, the gas abundance and the growth rate of the stellar 
component, 



eZ dZ s 
S„ dt 



(27) 



the constant e will, in the following, be treated as an essentially 
free (but not unconstrained) parameter of the model. The expected 
value is on the order of a few hundred, which is required to obtain 
r gr ~ 10 7 yr, suggested above. We will here adopt 



\ dt J \ Z) T T gr ' 



(28) 



as the rate of change of the dust-to-gas ratio Z d due to accretion 
of metals onto pre-existing dust grains in the ISM. Note that this 
fo rmulation of 'sec onda ry' dust product ion differs from that used 
by Edmunds j200 lh and lMattssonl ( 1201 lh in that it also depends on 
the dust abundance in the ISM and the depletion of metals in atomic 
state. 



efficiency is set by the fraction / d of interstellar dust destroyed in 
an encount er with a SN shock wave, which occurs in the definition 
of <mi SM > jMcKee 1989l: lDwek et~aT1l2007h . 



<miSM> 



•r 



/d(V.) 



dv„ 



dv s , 



(32) 



where M sw is the swept-up gas mass (during Sedov-Taylor ex- 
pansion), v s is the shock velocity, and vo, Vf are the initial and 
the final velocity, respectively. Note that in this way 6 is similar 
to the g-param eter (average grain-destruction efficiency) used by 
dMcKedfl989l) . which should not be confused with the e (dust- 
growth parameter) introduced in the previous section. It should also 
be stressed that f A is not a constant, but a function of the shock ve- 
locity iv 

A iLarsonl d 19981) IMF and <m ISM > ~ 1000M o dDwek et al.l 

l2007h suggests (5 ~ 10, which is likely close to an upper limit for 
6. Just as in the case of e above, it is not absolutely clear, however, 
that 6 can be treated as a parameter that does not vary during the 
course of evolution of the ISM in a galaxy, but it seems in a given 
environment a fair approximation. 



5 ANALYTIC SOLUTIONS 

For simplicity we have assumed a closed box (see section |3), i.e., 
no in- or outflows to/from the disc. This is not in agreement with 
the widely accepted ideas about galaxy-disc formation, where the 
baryons (in the form of essentially pristine gas) are assumed to 
be accreted ove r an extended period of time. But as shown by 
lEdmundsl dl990t) . the only major effect of unenriched infall is to 
make the effective yield smaller, i.e., to dilute the gas so that the 
metallicity builds up more slowly. As we in this study uses the 
present-day metallicity as input, the overall effects of assuming a 
closed box are rather small, and in general only accretion of metal- 
enriched gas can affect the dust-to-metals ratio significantly (see 
Appendix IBl. 



4.2 Dust destruction 

The dominant mechanism for dust destruction is by sputtering in 
the high-velocity interstellar shocks dri ven by SNe, which can b e 
directly rel ated to the ene r gy of the SNe l Nozawa & Ko zasa 2006). 
Following McKee (1989); iDwek et ail d2007l) the dust destruction 
time-scale is 



(mi SM )Rs 



(29) 



where £ g is the gas mass density, (otism) is the effective gas mass 
cleared of dust by each SN event, and R SN is the SN rate, which 
may be approximated as 



plOOAfo 



<p(m) dm. 



(30) 



The integral in equation yOJ is a constant with respect to time, and 
is not likely to vary much over the disc either, hence the time scale 
T d may be expressed as 

S dl. s 
\Z~dt' 



r- 1 « 



(3D 



where 6 will be referred to as the dust destruction parameter. This 
parameter is dimensionless, and as such it can be seen as a measure 
of the efficiency of dust destruction. More precisely, however, the 



5.1 General solution 

Adopting the closed-box scenario, the dust destruction and dust 
growth models as described above, results in an equation for dust 
evolution, 



* {* +2 #-t) z -']}lF' 

which combined with the metallicity Z gives 

i-sMK 1 -!)*-*!- 



(33) 



(34) 



where y z is the metal yield. Provided y A < y z , the general closed- 
box solution (of equation I34t for the dust-to-gas ratio Z d in terms 
of the metallicity Z is (see Appendix |C]f or a sketchy derivation), 



y±(z- 8 -\ 

yz V el 



yz 

for 

<Pij(k,x,y) = M 





y d (eZ-6) 2 s 2 ] 
yz' <yz ' tvz J 


-<pn 


.)' d 6 2 (eZ-6) 2 1 
yz ' eyz ' eyz \ 


2^io 


1" yj S 2 (eZ-6) 2 
[ yz ' eyz ' sj'Z 


+ Voi 


r.Vd (eZ-6) 2 S 2 1 

L yz ' tyz ' eyz J 



k 1 x 
i+ - , i + —; — 

2 2 2 



U 



k 1 y 
J 2 J 2 2 



(35) 



(36) 



The functions M and U are the confluent hypergeometric Kummer- 
Tricomi functions of the fir st and second kind, respectively 
dKummeill 1 837k iTricomUl 1 9471 . see also Appendix [Dj. 
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In the equations above, vj and y z are the stellar dust and 
metal yields, respectively, S is the 'dust destruction parameter' 
(see section |4~2t and e is the 'grain-growth parameter' (see sec- 
tion [47TJ. equation d35 b is singular at Z = 6/e, which means this 
general solution must be used with care. It is relatively straight 
forward to implement the Kummer-Tricomi functions numerically 
(see Appendix [Dj, but there is a regular singularity at the origin 
in Kummer's equation (to which M and U are linearly indepen- 
dent solutions) which can cause potential problems in the vicinity 
of Z = S/e. 



5.2 Special cases 

The general solution presented above is obviously not simple to 
use in practice, not the least because of the singularity at Z = S/e. 
However, in the special case y d — > (negligable net contribution 
of dust from stars) the singularity can be removed and the solution 
expressed as (see Appendix |Cl> 



-A = exp (£ 2 ) {exp + 77oZ d , [erfi<0 - erfi(f )]}~ 



with Zj o being the initial dust-to-metals ratio, 



(37) 



6 GRAPHIC ANALYSIS 

Using the numerical implementation of M and U described in Ap- 
pendix [D] we will here demonstrate the general behaviour of the 
dust-to-metals ratio £ = Z d /Z using contour plots. Unless anything 
else is stated, we assume y z = 0.02 is a good typical metal yield 
(which is consisten t with the results of Paper II, but note that, e.g., 
iGarnett et alj2002l finds a lower value) and y A = \ry z for simplicity. 



6.1 General effects of growth destruction of dust in the ISM 

In case of no dust destruction (6 = 0) the dust-to-metals ratio builds 
up to a maximum (where ( ~ 1) as e and the metallicity Z increases 
(see FigureQ] left panel). At low metallicities (half of solar or less, 
in the present case) the effect of increasing e is relatively small once 
we get beyond a certain e, while at higher metallicities f grows 
rapidly until the metals reservoir is exhau sted and ( approaches 
unity (as also foun d in the models by, e.g.,[Hirashita & Kuoi l201 il : 
lAsano et al]|2012h . In case of no dust growth (e = 0) the dust-to- 
metals ratio is on a steep "downhill slope" (approaching ( = 0) for 
essentially all metallicities and (5-values on the considered interval 
(see figure[2] left panel). Note that £ is very small at high metallicity 
if there is significant dust destruction. 



1 (eZ - S)- 
2 



ft 



1 (6Z - S) 2 

2 



'/o 



(38) 



and Z the initial metallicity. In the solution above, erfi(z) is the 
imaginary error function, related to the ordinary error function 
erf(z) as erfi(z) = — jerf(jz), where erf(z) is defined as 



erf(x) : 



\w Jo 



dt. 



(39) 



Physically, one may interpret this solution as describing the subse- 
quent evolution (where the dust contribution from stars may be con- 
sidered negligible) after an initial phase of metal enrichment and 
dust formation leading up to the point where Z = Zq and Z^Zd.o- 
It may not be entirely realistic, but it demonstrates the interstellar 
"battle" between growth and destruction of dust grains in a very 
nice way (see Section l63l and Figure(3](. 

Even when there is a significant net contribution from stars, 
we can still find simpler solutions for special cases. In case there is 
no dust destruction by SNe (6 = 0) the solution reduces to 



U(l + la3;^] 

2 _ >d V 2 >'z 2 2 \z ) ^ 

d "v z M(ia,l;i^ ' 
V 2 yz '2 '2 yz ) 



(40) 



and if there is dust destruction, but no grain growth in the ISM 
(e = 0), then 



6.2 Dust-to-metals gradients 

As shown by the theorem in section[3]the effect of dust destruction 
and dust growth on the dust-to-metals gradient in a galaxy disc is 
to make it steeper or flatter. The effect of growth and destruction of 
dust in the ISM can be illustrated in a more intuitive fashion if we 
consider the specific effects on a given metallicity profile. We here 
assume that metals in a disc follows an exponential distribution, 



Z(R) 



. R 

Zoexpj- — 



(43) 



where we set the central metallicity to Zq = 0.055 and the e-folding 
scale length R is set to be the unit for the galactocentric distance. 
The right panel of figure [TJ shows how dust growth creates a dust- 
to-metals gradient that falls of with galactocentric distance and be- 
comes increasingly steeper as e increases (for e = the gradient 
is flat). Similarly, the right panel of figure [2] shows how dust de- 
struction creates an inwards gradient, starting from a flat gradient 
for 6 = 0. 

In the context of dust-to-metals gradients as signs of either net 
dust growth or net dust destruction, one should as well note that if 
the metallicity gradient and the dust-to-gas gradients are essentially 
flat, it is more or less impossible to distinguish between pure stellar 
dust production (albeit with a high stellar dust yield) and scenario 
including dust growth and/or dust destruction in the ISM. 



ry yd 

Zd = — 



■ exp -5 — 

yz 



(41) 



If there is neither growth, nor destruction of dust in the ISM (e 
8 = 0), we have the trivial case 



yz 



(42) 



corresponding to pure stellar dust production and obviously a flat 
dust-to-metals gradient. All the special cases above evade the in- 
convenient singularity atZ = die. 



6.3 Dust growths vs. destruction 

Growth and destruction of dust must likely occur together. As we 
describe in Appendix Iclthe general solution with both e and 6 non- 
zero, has a singularity at Z = S/e, which makes the analysis of 
how growth and destruction compete somewhat complicated and 
not least limited. Hence, we will here consider the special (and not 
entirely realistic) case where stellar dust production is considered 
negligible (see equation [37} starting from a point in time when the 
metallicity Z = 1.0 ■ 10~ 5 and the dust-to-metals ratio £ = 0.5. In 
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figure [3] we show f as a function of e and 6 for a fixed present- 
day metallicity Z = 0.02. Increasing the efficiency of dust destruc- 
tion counteracts the dust growth, which is shown by the "downhill 
slope" towards high S and low e values. Clearly, a high efficiency 
of dust destruction is not likely if there is to be a significant net 
production of dust without invoking a ridiculously short time scale 
for the dust growth in the ISM. More precisely, it is required that 
eZ » 6, which in case Z = 0.02 and 6 = 10, would imply e » 500. 
With such a large e the typical growth-time scale is down to ~ 10 5 
yr or less. As we mentioned in section l4~Tl e should not exceed val- 
ues of a few hundred if the dust-growth time scale r gr is to be con- 
sistent with the suggested numbers for the local ISM of the Milky 
Way. It is quite possible that t e1 can be significantly shorter in, .e.g., 
a denser environment, but a deeper analysis of this goes beyond the 
scope of this paper. 



6.4 "Critical" metallicity for dust growth 

Just aslZhukovska. Gail & Trielofil fe008l) . lHirashita & Kud feOllh 
and lAsano et al.l J2012h . we find that there exist a "critical metal- 
liciy" Z crit where the dust-mass contribution from grain growth in- 
creases rapidly. But this rapid increase over orders of magnitude 
occurs only if the stellar dust yield y A is significantly lower than the 
metal yield (see figuref4] left upper panel). Moreover, Z crit depends 
somewhat on the dust-growth time scale (or e), which can be seen 
in figure[4](right upper panel). Hence, Z crit should not be viewed as 
a universal constant. In fact, a reasonable definition of Z crit would 
be the metallicity at which stellar dust production and the net dust 
growth in the ISM contribute equally to build-up of the interstellar 
dust component. In such a case, adopting the model used above, 



eZ, 6 



(44) 



If dust growth dominates over stellar dust production and dust de- 
struction in the ISM is negligible, i.e., if y d /e <k 1 and <5/e <k 1, 
then Zcrit ~ Zj, which suggest Z crit ~ y&lyz- At metallicities below 
this value, the dust evolution (as function of metallicity) should be 
essentially identical to the case of pure stellar dust production - 
without any growth or destruction of dust grains in the ISM. In the 
right upper panel of figure [4] Z clit = y^/yz is marked by a vertical 
dashed line. At lower metallicities all model curves are indeed the 
same. 

This critical metallicity Z clit has an interesting implication for 
dust-to-metals/gas gradients, as it predicts the existence of bends 
also in logarithmic slopes and the existence of a critical galacto- 
centric distance in between an inner and an outer "plateau" where 
£ is constant (see Fig[4] lower panels). This non-linear feature is the 
consequence of the interstellar dust-growth rate (see equation |28t 
having a non-linear (Z?) term. Although equation J24l > and equation 
d28t together represent one specific model, all models of interstel- 
lar dust growth will be non-linear as long as they depend on the 
amount of dust and "free" metals (not locked-up in dust) available. 



7 DISCUSSION AND CONCLUSIONS 

We have shown that dust destruction by shock waves from explod- 
ing SNe and interstellar dust growth acts in opposite ways on the 
dust-to-metals gradient over a galaxy disc (see the theorem proved 
in section[3}- This is hardly surprising, but starting from an exactly 
flat gradient (or no gradient, more precisely) dust destruction will 
over time create an inwards slope, while dust growth will create an 



outwards slope, provided the dust-to-gas ratio as well as the metal- 
licity have negative gradients, i.e. decreases with galactocentric dis- 
tance. Hence, we expect dust-to-metals gradients to have positive 
(inwards) gradients if dust destruction is more important than dust 
growth, and if dust growth is the more important process we expect 
them to be negative in general. The dust-to-metals gradient thus ap- 
pears to be a useful diagnostic for the existence of interstellar dust 
growth. 

Our simple model of dust growth has just one adjustable pa- 
rameter. This parameter (e) can have a rather wide range of numeri- 
cal values depending on what one assumes about the physical prop- 
erties of the dust grains as well as the gas in ISM. In principle e is 
proportional to the gas mass density E g if all other quantities remain 
constant, but the star for mation efficiency n is li kely proportional to 
E e raised to some power (Krumholz & McKee 2005) and since the 
cloud size d c is also likely related to E g , e may not be much depen- 
dent on X g after all. More precisely, the star-formation efficiency (or 
time scale) is expected to correlat e with the free-fall time scale, i.e., 
E» oc I.fy/Tfi, where Tg oc E^ 1/2 jKrumholz. McKee & Tumlinsor] 
2009), assuming E g oc p B telmegreerj|2002l) . and the scale of the 
cloud size d c is given by Jeans length A s oc E g . As e oc t s d c E e 
and (v g ) is roughly constant (isothermal conditions), the effective 
dependence on E a is expected to be weak, if not negligible. Thus, 
it is fair to assume that e is (effectively) only very weakly depen- 
dent on Eg within a galaxy, although from one galaxy to another e 
may vary significantly, however (see Paper II). Below we analyse 
the range of possible e values considering just mean/characteristic 
values of Eg, ^ and <v g ). 

In terms of the included physical parameters (see section l4~TT l. 
we find 



f,n<r(v g ) 

£ » (Eg). 

aim («c> 



(45) 



The lock-up fractio n a is 0.6 - 0. 8 for a normal IMF (we use 
here a = 0.7, see lMattssonll201l[ figure 1), (77) is ~ 1 Gyr 1 
and since the typical size of a molecular cloud d c is 10 - 100 pc, 
we adopt (d c ) = 50 pc. The average grain mass <m gr ) of course 
depends on the typical grain size a, where the latter ranges be- 
tween 0.001/um for the smallest seed particles and ~ lyum for large 
full-grown dust grains. Hence, it is more convenient to introduce 
the characteristic grain density p gr = (m gI )/{V a ), where V a - is 
the vol ume of a dust grain . The grain density p gl is typically 3.3 
g cm 3 jDraine & Lill2007f) for silicates and 1.85 g cnT 3 for amor- 
phous carbon jRouleau & Martir]|l99ll) . but other values can also 
be found in the literature. Taking p gr = 2.5 g cnr 3 as representative 
figure for cosmic dust in general, we arrive at 



' /s Um/ \MePC- 2 



(46) 



Assuming that all metals that come in contact with a dust grain 
will stick to that dust grain (f s = 1), a small characteristic grain 
size a - O.Olyum and a relatively high average gas density of 
E 5 = 50 M G pc~ 2 , will result in an e of roughly 2- 10 4 corresponding 
to a typical grain-growth time scale of r gr ~ 10 6 yr if the gas con- 
sumption rate is similar to that of the solar neighbourhood. Such 
high values of e may be expected in young star-forming systems 
(e.g., late-type dwarf galaxies) where one has reasons to believe 
that gas densities are quite high and the grain-size distribution is bi- 
ased towards small grains. The latter is due to insufficient time for 
extensive grain growth, and grain shattering owing to an elevated 
SN rate and strong UV radiation as consequences of recent star for- 
mation. If f s = 0.1 (which is more consistent with silicate growth), 
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Figure 1. Left: Dust-to-metals ratio £ = Zj/Z as a function of the metallicity Z and the dust-growth parameter e for the case where there is no dust destruction 
due to SNe (<S = 0). Right: Same as the left panel, but as a function of the galactocentric distance in a galaxy disc assuming an exponential distribution of 
metals. 
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Figure 2. Left: Dust-to-metals ratio £ = Z&/Z as a function of the metallicity Z and the dust-destruction parameter S for the case where there is no dust growth 
in the ISM (e = 0). Right: Same as the left panel, but as a function of the galactocentric distance in a galaxy disc assuming an exponential distribution of 
metals. 



a = lyum and E e = 5 M e per - , then e ~ 2 and thus t si ~ 10 yr. 
With f s = I, a = 0. lyum and E e = 10 M e per 2 , the corresponding r sr 
is ~ 10 8 yr assuming a gas-consumption rate and a metal content 
similar to that of the local Galaxy. Thi s number is consistent with, 
estimates made in some other stud ies (Jones, Tielens & Hollenbach 
1996t |jonesll2004 IZhukovska. Gail & Trielofll2008M Jones & Nutb! 



201 1[), but slightly longer than the time scales suggested recently 
bv lHirashita&Kud feOllh . Note that grain size is the parameter 
that is likely most important for the value of e as it can vary quite 
significantly. The gas mass density can vary over several orders of 
magnitude as well, but as described above it may be cancelled out 
by other parameters. 



The model of dust destruction due to SN shock waves has ef- 
fectively only one parameter as well. This dust destruction param- 
eter 6 can be expressed as 



r* J UU/l 



(p(m) dm, 



(47) 



where cf> is the IMF and a is the lock-up fraction, as previousl y 
defined. With a = 0.7 and a normal IMF (see, e.g., lLarsonlll998r) . 
we find 



(48) 



which suggest 6 is of order ten, if (/7Jism) ~ 1000 M Q . The actual 
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Figure 3. Contour plot of the dust-to-metals ratio as function of the parameters e (growth) and S (destruction) for the special (and not entirely realistic) case 
where stellar dust production is considered negligible (see equation |37l starting from a point in time when the metallicity Z = 1.0 ■ 10~ 5 and the dust-to-metals 
ratio I = 0.5. The present-day metallicity is assumed to be roughly solar (Z = 0.02). 



efficiency of dust destruction, and thus the effective interstellar gas 
mass cleared of dust, is not very well known. Therefore, it is reason- 
able to treat 6 as an essentially free parameter. In order to have net 
growth of dust in the ISM, the value of 6 needs to be 6 < e (Z — Zd). 
This means S ~ 10 is likely at the upper end of possible values for 
such a scenario, assuming Z - « 0.01, which suggest 6 = 10 
would require e > 1000 or r gr < 10 7 yr. High values of d may be 
found in starburst environments, where high SN rates and possibly 
also top-heavy IMFs are expected. However, in general 5 is likely 
small, since high rates of dust destruction are somewhat inconsis- 
tent with the fact that dust is ubiquitous throughout the Universe. 

Although simplifying assumptions have been made in this 
study in order to obtain a reasonably simple parametric model in 
terms of e and 6, a clear outwards slope is unlikely to be the re- 
sult of any other mechanism than dust growth in the ISM. Other 
mechanisms, which however appear less effective: 

• Accretion of dust free material onto the galactic disc may 
affect the dust-to-metals ratio if the infalling gas contains some 
fraction of atomic metals (see Appendix |B1 for further details and 
worked out examples). The metallicity of the accreted gas is likely 
much less than that of the ISM, so the effect cannot be very large 
and it would also mimic the effect of dust destruction rather than 
dust growth. 

• Secondary dust production in stars, i.e., a stellar dust yield 
which increases as the metallicity of stars increases, may in princi- 
ple create a dust-to-metals gradient along a galaxy disc. However, 
the relative increase of the stellar dust yield along the disc cannot 



be arbitrarily large. In particular, the dust-to-metals gradient can 
never become steeper than the metallicity gradient only owing to 
secondary dust production in stars (for further details and a more 
quantitative analysis, see Appendix lAl. 

• The lifetime of stars may also play a role, but since the very 
same stars that are producing the metals are also responsible for the 
stellar production of dust, this effect cannot be dominant. In fact, it 
should be negligible. 

Thus, we conclude that dust-to-metals gradients can be used 
as a diagnostic for interstellar dust growth in galaxy discs, where 
a negative slope indicates dust growth. 

Dust growth has a non-linear nature as the time scale for 
it must depend on both the metallicity and the amount of avail- 
able seed grains. As a consequence there is a "critical" metallic- 
ity (which depends on the dust-growth and dust-destruction time 
scales as well as the dust-to-gas ratio) at which the dust production 
by interstellar grain growth exceeds stellar dust production and the 
dust-to-gas ratio diverges from the steady increase obtained in case 
the dust mass is owing to stars only. This allows for bends in the 
logarithmic slopes of the dust-to-metals profile even if the metal- 
licity follows an exponential fall-off with galactocentric distance. 
Dust destruction in the ISM due to SNe may also affect the shape 
of the dust-to-metals profile, creating a central depression as the 
dust-to-gas ratio, the metallicity and the integrated number of SNe 
typically increases in the central parts of a galaxy disc compared 
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Figure 4. Effects of the critical metallicity for dust growth domination and its dependence on the stellar dust yield and dust growth parameter e. The upper 
panels show the evolution of the dust mass as a function of metallicity for various values of the stellar dust yield y$ with a fixed e = 200 (left panel) and 
various values of e with a fixed stellar dust yield yj = 5.0 ■ 10~ 3 yz (right panel). The lower panels show the corresponding plots of the dust-to-metals ratio ( 
as a function of galactocentric distance assuming an e-folding decay of the metallicity along the disc (see equation pBV 



to the outer disc. However, since dust growth increases as well, the 
expected net effect is an increased dust-to-metals ratio in any case. 

Finally, we note that combining recent obser vational results 
jMunoz-Mateos et alj|2009l ; iMoustakas et al.ll201oh one finds that 
dust-to-metals gradients in late-type galaxy discs appear relatively 
steep (and negative), i.e., show a clear fall-off with galactocentric 
distance, which suggest interstellar dust growth is more important 
than stellar dust production. In Paper II of this series, where we 
compare theoretical models and observational results in more de- 
tail, we return to this fact and look for more quantitative evidence of 



interstellar dust growth being the dominant dust production mech- 
anism in late-type galaxies. 
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APPENDIX A: EFFECTS OF SECONDARY STELLAR 
DUST PRODUCTION 

It is important to remember that interstellar dust growth is not the 
only mechanism that can give rise to a dust-to-metals gradient. Sec- 
ondary dust production in stars, i.e., a metallicity-dependent yield, 
could in principle have similar effects. Splitting the stellar dust 
yield into two components, the constant primary yield yj! and the 
metallicity-dependent secondary yield y s A = yl Q Z/Z Q , and assum- 
ing there is no growth, nor destruction of dust in the ISM, we obtain 



dZ 



1 

yz 



(Al) 



which, with the initial condition Zd(0) = 0, has the solution 
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z / p y d ,o z 



yz V " 2 Z c 
Using the notation introduced in section[3] we can also write 



4 I -^.0 Z 



(A2) 



(A3) 



In the limit where y|j <sc yj Q Z we then have ~ A z , which is the 
steepest dust-to-metals gradient obtainable for a given metallicity 
gradient. Hence, if the dust-to-metals gradient is steeper than the 
metallicity gradient, then there must be dust growth in the ISM to 
account for that steepness. 

It is quite unlikely that yjj <K yjj Z as dust production in stars 
is likely primary to almost the same extent as the metals production 
is. Since most of the metals are primary, the secondary dust compo- 
nent cannot be dominant, which implies A ( < A z . In fact, modelling 
of stellar dust production suggest dust yields can b e relatively high 
even at Z = (see lGal l. Hiorth & Anderser] |201ll . and references 
therein). It is actually reasonable to assume the secondary yield is 
no more (likely less) than 50% of the primary yield at solar metal- 
licity Zq. If y^ Q 



fj% then 



(A4) 



where we note that 4/Z G ~ 300. More precisely, this implies A( is 
at least about two orders of magnitude smaller than A z for all rea- 
sonable metallicities along a galaxy disc. Thus, we conclude that al- 
though the dust-to-metals gradient A f is technically non-zero in this 
case, it is still consistent with a flat du st-to-metals profile as metal- 
licity gradients are rarely very steep jPilvugin. Vflchez & ContinU 
2004). Secondary dust production in stars cannot be responsible 
for a significant dust-to-metals gradient. 



APPENDIX B: EFFECTS OF INFALL 

Throughout this paper we have treated the dust evolution in late- 
type galaxies assuming they are "closed boxes", i.e., that there is 
neither any inflow, nor any outflow of gas and metals to/from the 
disc. In reality accretion of gas and minor mergers with smaller 
galaxies are important for the chemical evolution of a galaxy disc 
and thus also important for the shaping of the dust component. 
Hence, an infall component in equation < |34t would have been in its 
place, but we omitted it for simplicity. However, as we will show 
here, the effect of infall is not such that it can qualitatively change 
any of our results. In fact, a theorem similar to that presented in 
section [3] could likely be formulated, but it would be less transpar- 
ent as regarding the effects of growth and destruction of dust in the 
ISM. 

In case of no growth or destruction of dust in the ISM and 
accretion of pristine (unenriched) gas at a rate % n {, which contains 
no metals in any form, we have the equation 



dZi 

dZ y z - ZA 



(Bl) 



If A is constant, the solution to the equation above is 

Z d = — Z, (B2) 

yz 

given the initial condition Zd(0) = Z(0) = 0. Hence, pristine infall 
likely does not affect the dust-to-metals ratio much. 

If the accreted gas contains metals the situation is quite differ- 
ent. Including metal-enriched infall, equation iB It becomes 



dZ A 

Hz 



Z,A 



(B3) 



yz - (Z-Z inf )A 

Assuming an outflow of interstellar gas where some of the metals 
in that gas are accreted back onto the disc, we can wr ite Z m f = vZ 
(usually referred to as a "galactic fountain" model, see lRecchi et al.1 
2008) the solution is 



1 + (v- 1)A — 



(B4) 



where y' z is a reduced metal yield to account for the metals lost in 
the outflow. For the specific case v = 1/2, we obtain the solution 



Zi 

z 



Id 

y 7 



i - 



A Z 

4yT 



(B5) 



from which it is easy to see that the effect of infall is reminiscent 
of the effect of dust destruction in the ISM due to SNe. i.e., that 
the dust-to-metals ratio decreases with metallicity (cf. equationl4"Tll. 
Moreover, if the dust destruction term is included in equation IB lb . 
we have 



yd 

A + S 



1 + (y- 



»4 

y z 



(B6) 



which is a solution of the same mathematical form as above. Hence, 
it is quite clear that metal-enriched infall has an effect which is very 
similar to that of dust destruction by SNe, which means that infall 
alone cannot create a dust-to-metals gradient with the same sign as 
the metallicity gradient. 



APPENDIX C: GENERAL SOLUTION OF EQUATION (34) 

The general solution to equation d34t presented in section [5] is 
expressed in terms of a product of the confluent hypergeomet- 
ric K ummer-Tricomi functions M and U jKumme3l 1 8371 : iTricomil 
1 1947b . This solution exist because equation d34t is related to Rum- 
mer's equation (also known as the confluent hypergeometric equa- 
tion), i.e., 

d 2 w dw 
z— T +(b-z)—-aw = 0, (CI) 

dz dz 

which has the solution w{z) = c\M(a,b;z) + C2U(a,b;z), where C\ 
and c'2 are arbitrary constants. With the variable change 

1 (eZ - S) 2 

-. C2) 

2 ey z 

equation d34b can be rewritten as an equation of the form 
dZ A _ y d „ e 



+ Z d 



(C3) 



This is a Riccati equation, which has the general form 

% = q (x) + c h {x)y(x) + q 2 (x)y 2 (x). (C4) 
dx 

Such non-linear equations ca n be reduce d to a second order linear 
ordinary differential equation dlncell 19561) of the form 

d 2 u du 

— - R(x) — +S(x) u(x) = 0, (C5) 
dx 1 dx 

where 

R(X) = qi(x) + -^—^, S(x) = q (x)q 2 (x). (C6) 
q2\x) dx 

A solution to equation JC5b provides a solution to equation l lC4b as 
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y{x) = — 



l 



du 



q2(x) u(x) dx 
Identifying the Riccati coefficients as 

yd 



9>(g) : 



h(f) = l, 92(f) = — 



(C7) 



(C8) 



we find the associated second order linear ordinary differential 
equation, 

d 2 u 1 1 \ 1 y,i 

\ 2 / d£ 2 y z 

This is the Kummer equation for b = 1/2 and a = y^/2y z . Reverse 
Riccati reduction and back-substitution, together with the natural 
initial condition Zd(0) = 0, will provide the general solution given 
in section[5]after some algebra. 

The Kummer equation has an awkward property: it has a regu- 
lar (order one) singularity at the origin (at f = in the case above). 
This means that no solution exists at £ = (or Z = 6/e) and that 
the region near this point must be avoided when applying this so- 
lution. In particular, when the Kummer- Tricomi functions are im- 
plemented numerically, the algorithm for computing them will be 
unstable in the vicinity of this singular point. We have therefore 
considered the special (and not entirely realistic) case where stellar 
dust production is negligible from a point in time when the metal- 
licity has a certain value Z and the dust-to-metals ratio has a value 
fo. Assuming y<j — > equation Q4t reduces to 



yz 



dZ A 
dZ 



.(eZ-S)Z i -eZ\. 



Riccati reduction as above yields 



du 

-rr =0, 



(CIO) 



(Cll) 



2-t) d$ 

where f is as previously defined (equation IC2t . This equation is 
non-singular, but requires that Z > Z , where Z is some finite ini- 
tial value. The general solution is 



«(f) = C + C 1 V5rerfi(VI), 



(CI 2) 



where Co, Ci are constants to be fixed by initial conditions as we 
do reverse the Riccati reduction and back-substitute. With Zo = 
Z(f ) t and Z d = Z d (0) ^ as the initial conditions it is then 
possible to obtain the solution given as equation i37l in section 
(ED. 



where 



: a(a + l)(a + 2) • • • (a + n — 1) 



T{a + ri) 



T(a) ' 

is the Pochhammer symbol and T is the Gamma function, 



r(z) = f / - l c~'dt. 

Jo 

The function M can thus be implemented numerically by comput- 
ing the above series until some arbitrary precision is obtained. The 
typical number of terms needed to reach the precision limit of a 
standard Intel processor is at most a few hundred. This may still 
cause problems when computing the Pochhammer symbol, since 
this will have to be done using some limited implementation of T 
to obtain reasonable computation speed. The basic issue is the fact 
that T, as well as the factorial, is usually not implemented for large 
arguments. For example, in IDL and MATLAB the argument z can- 
not exceed ~ 170. However, this situation rarely occurs. 

The func tion U can be denned in terms of the function M 
jTricomill 19471) by 

r(i - b) 



(D3) 



(D4) 



(D5) 



U(a,b;z) : 



T(a - b + 1) 



M(a, b; z) + 



(D6) 



^-Jl z ^M(a-b+ 1,2 -b;z), 
T{a) 



which is straight forward to implement, except for integer b (where 
U is not defined). The general solution to equation J34t does not 
give rise to any integer values for the b, so we will not consider 
how to implement the analytical extension of (J for integer b. This 
can be done, but goes beyond the scope of this study. 



APPENDIX D: NUMERICAL IMPLEMENTATION OF 
THE KUMMER-TRICOMI FUNCTIONS 

The Kummer- Tricomi functions, used in section [5] and Appendix 
|Cl above, can be defi ned in terms of integral quantities jKummeJ 
[l837tlTricornl 1947b . 

M(a,b;z) = —^ - f e z "u a -\i-u) b - a - l du, (Dl) 
T(a)T(b - a) J 

for %(b) > %(a) > 0, 

U(a,b;z)= — f e n f^{\ + t)"'"' 1 dt, (D2) 
r(a) Jo 

for Dl(a) > 0. This is not very convenient for numerical implemen- 
tation though. As alternative we can consider the following. The 
function M is identical to the iFj hypergeometric function which 
can be defined as an infinite series, 
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